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We use homological perturbation machinery specific for the algebra category 
1 13] to give an algorithm for computing the differential structure of a small 1- 
homological model for commutative differential graded algebras (briefly, CDGAs). 
The complexity of the procedure is studied and a computer package in Mathematica 
is described for determining such models. 



1 Introduction 



The description of efficient algorithms for homological computation can be considered 
to be a very important topic in Homological Algebra. These algorithms can be used 
mainly in the resolution of problems in Algebraic Topology; but this subject also impinges 
directly on the development of diverse areas such as Combinatorial Designs, Code Theory, 
Concurrency Theory or Cohomological Physics. 



Starting from a finite CDGA A, we establish an algorithm for obtaining an "eco- 
nomical" 1-homological model hBA, in the sense that the number of algebra generators 
of hBA is less than that of the reduced bar construction B{A). In order to get the 
1-homology of A, we would need to compute the homology groups of the model hBA. 

*Proceeciings of the 3rd Workshop on Computer Algebra and Scientific Computing CASC'OO, 
Springer-Verlag (2000), 87-100. 

'I' Alvarez V., Armario J. A., Frau M.D., Gonzalez-Diaz R., Jimenez M.J., Real P. and Silva B. Authors 
are partially supported by the PAICYT research project FQM-0143 from Junta de Andalucia and the 
DGES-SEUID research projects PB98-1621-C02-02 from Education and Science Ministry (Spain). 
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This computation can be reduced to a simple problem of Linear Algebra (see |TT|] for a 
complete explanation of this method). 



Our main technique is homological perturbation machinery ^. Homological 

Perturbation Theory is often used to replace given chain complexes by other smaller, 
homotopic chain complexes which are more readily computable. An essential notion in 
this theory is that of contraction. A contraction c = {f,g,(j)) between two differential 
graded modules (iV, cIm) and (M, (Im) is a special homotopy equivalence between both 
modules such that the corresponding homology groups are isomorphic. The morphisms 
/, g and are called projection, inclusion and homotopy of the contraction, respectively. 
The Basic Perturbation Lemma is the heart of this theory and states that given a con- 
traction c = {f,g,(f)), and a perturbation 6 of d^^ (that is, {dj^ + 6)"^ = 0), then there 
exists a new contraction cs = {fs, gs, 'Ps) from (A^, d^^ + 6) to (M, d^i + ds), satisfying 

/, = /(I - , gs = K9, 05 = Sf0 

ds = f5J:lg, (1) 
where = ^(-1)^(05)^ = 1-05 + 0505 + (-1)^(05)* + ■■ -. 

j>0 

It is necessary to emphasize that a nilpotent condition for the composition 50 is 
required for guaranteeing the finiteness of the formulas. 

The basic idea we use in this paper is the establishment (via composition, tensor 
product or perturbation of contractions) of an explicit contraction from an initial differ- 
ential graded module A^ to a free differential graded module M of finite type, so that the 
homology of A^ is computable from that of M. 

This "modus operandi" has been used by the authors in previous works [0, |, 3. 

Working in the context of CDGAs, Homological Perturbation Theory immediately 
supplies a general algorithm computing the 1-homology of these objects at graded mod- 
ule level. Nevertheless, this procedure, already presented by Lambe in [|r^, bears, in 
general, high computational charges and actually restricts its application to the low 
dimensional homological calculus. 

This algorithm is refined, taking advantage of the multiplicative structures, in p|. 
More precisely, the Semifull Algebra Perturbation Lemma |T3|, Sec. 4] is used for design- 
ing the algorithm Algl. The input of this method is a CDGA A given in the form of 
a "twisted" tensor product of n exterior and polynomial algebras, and the output is a 
contraction cs (produced via perturbation) from the reduced bar construction B{A) to a 
smaller differential graded algebra hBA, which is free and of finite type as a graded mod- 
ule. In this case, we say that the pair {cs, hBA} (or, simply, hBA) is a 1-homological 
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model of A. Taking advantage of the fact that the differential d^BA of hBA is a deriva- 
tion, that is, a morphism compatible with the product of the 1-homological model, it 
is only necessary to know the value of this morphism applied to the generators of the 
model (let us observe that there are n algebra generators). This implies a substantial 
improvement in the computation of the differential on the small model hBA. 

In spite of this improvement, the computational cost for determining the morphism 
dhBA applied to an algebra generator of hBA is enormous, since the differential d^BA 
follows the formula (|l]) and the homotopy of c has an essentially exponential nature 
not only in time but also in space. 

We develop some techniques, which comprise what we call Inversion Theory and 
which first appears in |[T^. In consequence, we refine the formula for which is involved 
in the description of d^BA- This study is based on the observation that the the projection 
/ applied to certain elements (those "with inversions") is always null. It follows that a 
not insignificant number of terms in the formula of the morphism cj) can be eliminated 
in the composition f6{(j)6yg, which appears in the formula (|l]) of ds. In such a way, we 
derive an upgraded algorithm Alg2. 

The article is organized as follows: Notation and terminology are introduced in Sec- 
tion ^. In Section ^ the algorithm Algl, which was described in is recalled. Our 
contribution starts in Section ^ which is devoted to explaining Inversion Theory and 
describing the algorithm Alg2. An analysis of the complexity of Algl and Alg2 for 
computing the differential structure of the small 1-homological model hBA is carried 
out in Section ^ and a comparison between both algorithms is given. Finally, in Section 
^ we also give several examples illustrating the implementation of Alg2 carried out using 
Mathematica 3.0. 



2 Preliminaries 



Although relevant notions of Homological Algebra are explained through the exposition 
of this paper, most common concepts are not explicitly given (they can be found, for 
instance, in [|10] or |ll5|). 



Let A be a commutative ring with the non zero unit, which will be considered to be 
the ground ring. A DGA-module {M,dni,^M,rin.j) is a module endowed with: 

• A graduation, that is, M = ©„gN^n- 

• A differential, du '■ M ^ M, which decreases the degree by one and satisfies 
dl = 0. 

• An augmentation, : M — >• A, with ^m<^i = 0. 
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• A coaugmentation, rj : A —>■ M, with C,MriM = 1 



We will respect Koszul conventions. The homology of a differential graded module M, is 
a graded module H^{M), where Hn{M) =KeT dn/lmdn+i- We are specially interested on 
CDGAs, {A, d^, *A, ^A, Va) which are differential graded modules endowed with a product, 
*A, that is commutative in a graded sense. A morphism 6 : A ^ A which decreases the 
degree by one, is a derivation if 6* a = (g) 5 + 5 (g) 1). 

Three particular algebras are of special interest in the development of this paper: 
exterior, polynomial and divided power algebras. Let n be a fixed non-negative integer. 



• The exterior algebra E{x, 2n + 1) is the graded algebra with generators 1 and x of 
degrees and 2n + 1, respectively, and the trivial product, that is, x • x = and 
X ■ 1 = X. 

• The polynomial algebra P{y, 2n) consists in the graded algebra with generators 1 
of degree and y of degree 2n. The product is the usual one in polynomials, i.e.: 
y^ ■ y^ = y*'^"', for non negative integers i and j. 

• Finally, the divided power algebra T{y, 2n) is the graded algebra with generators 1 
and y { y = y^^^) of respective degrees and 2n. The product is defined by the 

rules T/^*-* ■ y'^^^ = i ^ ^, j y^'^^^\ i and j being non-negative integers. 



Each one of these three types of algebras can be considered as a CDGA with the trivial 
differential. 



Now, we shall recall a standard algebraic tool which allows us to preserve the product 
structure of the initial CDGA through the procedure of homological computation. The 
reduced bar construction |T0[ associated to a CDGA A is defined as the differential graded 
module B{A): 



B{A) = A © Ker © (Ker © Ker © ■ ■ ■ © (Ker G ® ■ ■ • ® Ker © ■ ■ ■ . 

An element from B{A) is denoted by a = [ai| ■ ■ ■ |a„]. There is a tensor graduation 
I \t given by ■ ■ • |cin]|t = YJi=i kil; as well as a simplicial graduation \ |s, which is 
defined by \a\s = |[ai| ■ ■ ■ |a„] 1^ = n. The total degree of d is given by \d\ = \a\t + \a\s. 

The total differential is given by the sum of the tensor and simplicial differentials. 
The tensor differential is defined by: 

dt[ai\ ■ ■ ■ \an] = - X:(-l)l'"^l-l"'-^ll[ai| ■ ■ ■ M^a,| ■ ■ ■ |a„] . 
1=1 
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The simplicial differential acts by cutting down the simphcial degree by using the product 
given in A. 



When the algebra A is commutative, it is possible to define a multiplicative structure 
on 13 (A) (via an operator called the shuffle product), so that the reduced bar construction 
also becomes a CDGA. 

Given two non-negative integers p and q , a {p, q) -shuffle is defined as a permutation 
TT of the set {0, . . . ,p + g — 1}, such that 7r(i) < 7r(j) when Q<i<i<p — lox 
p<i<j<p + q~l is the case. 



Let us observe that there are 




different (p, g)-shuffles. 



So, given a CDGA A, the shuffle product i< : B{A) (g) B{A) — > B{A), is defined by: 
[ai\ ■ ■ ■ \ap\ -k [bi\ ■ ■ ■ \bg] = J2 (-l)^^^'"'*^[c7r(o)|---|c^(p+g-i)]; 

Te{(p,i3) shuffles} 

where (cq, . . . , Cp_i, Cp, . . . , Cp+g_i) = (ai, . . . , Op, 6i, . . . 6,) and 

7r(i)>7r(p+j) 

Let be a non-negative integer. The ?7,-homology of a CDGA A (see |T^) consists 
in the homology groups of the iterated reduced bar construction B"-{A) = B{B"^^{A)), 
being B°{A) = A. 

Let {Ai}i^j be a set of CDGAs. A twisted tensor product ^^^jAi is a CDGA satisfying 
the following conditions: 

i) ^i^jAi coincides with the tensor product ^i^jAi as a graded algebra. 

ii) The differential operator consists in the sum of the differential of the banal tensor 
product and a derivation p. 

A contraction c : {N, M, f, g, (j)} |^, also denoted by {f,g,(p) : 4> M, from 
a differential graded module (A^, d^f) to a differential graded module (M, rf^/) consists 
in a homotopy equivalence determined by three morphisms /, g and f : ^ 
(projection) and g : M^, — > A^^, (inclusion) being two differential graded module 
morphisms and (j) : A^*+i a homotopy operator. Moreover, these data are required 

to satisfy the following rules: 

fg = 1„ , (pdj^ + d^cj) + gf = 1^, /</> = , (j)g = 0, </)</) = . 
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There are two basic operations between contractions which give place to new con- 
tractions: tensor product and composition of contractions. 



In this paper we use a particular type of contraction between CDGAs. Given two 
CDGAs A and A', a semifull algebra contraction {f,g,(f)) : A ^ A' consists of an 
inclusion g that is a morphism of CDGAs, a quasi-algebra projection / and a quasi- 
algebra homotopy (p. We recall that 

1. The projection / is said to be a quasi-algebra projection whenever the following 
conditions hold: 

/(0*^0) = O, /(0*^^) = O, f{g*^<l)) = 0. (2) 

2. The homotopy operator (f) is said to be a quasi-algebra homotopy if 

^0) = O, 0(0*^^?) = O, cp{g*^<p)=0. (3) 



The class of all semifull algebra contractions is closed under composition and tensor 
product of contractions. Moreover, this class is closed under perturbation. 



Theorem 2.1 

Let c : {N, M, f, g,(j)} be a semifull algebra contraction and 6 : N ^ N be a 
perturbation-derivation of d^- Then, the perturbed contraction cs, is a new semifull 
algebra contraction. 



3 Computability of the 1— Homology of CDGAs. 
First Algorithm 

Here we recall the algorithm described in for the computation of a 1-homological 
model of a CDGA. 

It is commonly known that every CDGA A "factors", up to homotopy equivalence, 
into a tensor product of exterior and polynomial algebras endowed with a differential- 
derivation; in the sense that there exists a homomorphism connecting both structures, 
which induces an isomorphism in homology. 

In fact, our input is a twisted tensor product of algebras A = ^'^^jAi where / denotes 
a finite set of indices, p is a differential-derivation and Ai an exterior or a polynomial 
algebra, for every i. In our algorithmic approach, we encode A by 
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1. a sequence of non-negative integers rii < n2 < ■ ■ ■ < Uk, such that rii represents 
the degree of the algebra generator Xi of Af, 

2. a /c-vector v = {vi, t>2, . . . , Vk), such that Vi is p(xj) for all i. 

The principal goal is to obtain a "chain" of semifull algebra contractions starting at 
the reduced bar construction B{A) and ending up at a smaller free (as a module) CDGA. 
In that way, we determine a 1-homological model for A. 

Now, we consider the following three semifull algebra contractions which are used, 
firstly, to find the structure of a graded module of a 1-homological model for a CDGA: 

• The contraction defined in [g § from B{A ® A') to B{A) ® B{A'), where A and 
A' are two CDGAs. 

Cs« : {B{A ® A'), B{A) ® B{A'), f^^, g^^, 0^^} ; 

- /s®[ai ® a'll ■ ■ ■ |a„ (g) a^] 

n 

= XI ^-4(0,+! *A ■ ■ ■a„)G'(«'i *A' ■ ■ ■a')[ai| ■ ■ ■ |aj O [a'^J ■ ■ ■ \a'J\ 

- ■ ■ ■ kn] ® Kl ■ ■ ■ laj) 

= [ai ®e'\--- |a„ ® * ® a'J ■ ■ ■ 1^ ® a^] , 

where 6 and ^' are the units in A and A' respectively. 

- up to sign, (j)B^{[ai a'^l ■ ■ ■ |a„ (g) a'J) 

= ±G(«n-,+i *A ■ ■ ■a„)[ai a'l ■ ■ ■ |a^_i (g) a^_, 

where n = n - p - q, (co, . . . , c^+J = (a^, . . . , a„_„ ---a'J and the sum 

is taken over all the {p + 1, g)-shuffles n and 0<p<n — q — l<n — 1. 

Let us note that the complexity of gs^ and 0^® is exponential since shuffles 
are involved in both formulas. 

Given a tensor product (gjg/Aj of CDGAs, a contraction from i?((g)jg/y4j) to 
(g)jg/i?(y4j) is easily determined by applying C^^ several times in a suitable way. 
This new contraction is also denoted by C^®. 



7 



• The isomorphism of differential graded algebras (therefore, a contraction) 

CsE ■■ {B{E{u, 2n + 1)), r(M, 2n + 2), /^^, g^^, 0} 
described in where 

. ,r I m times , (^\ , (m\\ r i "i times , -, 

■■■ |m]) = m^ ^ gBE{u^ ') = [u\ ■■■ \u]. 

• The contraction 

C^P ■■ {B{P{v, 2n)), E{v, 2n + 1), /^p, Qsp, 0bp} 
stated in 0, where 

gBp{v) = [v] and 0sp([t;''n ■ ■ ■ l^^"""]) = [v\v''^''\ ■ ■ ■ . 

Thanks to these three contractions, it is possible to establish, by composition and 
tensor product of contractions, the following semifuU algebra contraction C = {f,g,(j)): 

where hBAi represents an exterior or a divided power algebra with a generator Xj, 
depending on whether Ai is a polynomial or an exterior algebra with a generator Xi. 

In order to obtain the differential structure of the 1-homological model for the twisted 
tensor product ®1^jAi, the next step is to perturb C . The perturbation p produces a 
perturbation-derivation 5 on the tensor differential of B{®i^iA.;): 

n 

•^IKI ■ ■ ■ W\) = ■ ■ ■ lp(«OI ■ ■ ■ K] . 

1=1 



Now, by applying Theorem |2.1| , a new semifuU algebra contraction {fs, gs, 4>s) is 
constructed: 



where the differential ds is determined by the perturbation procedure (Basic Perturbation 
Lemma). That means that hBA = ^i^j{hBAi,ds) is a 1— homological model of A = 
®jg/y4j. Let us emphasize that the Basic Perturbation Lemma provides finite formulas. 
Indeed, this is a consequence of two facts: the perturbation 6 does not change the 
simplicial degree and increases this degree. 
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Procedure 1 Algorithm Algl. 



Input; A finite CD G A A: {{rii, . . . ,12^), {vi, . . . ,Vk)) ■ 

Output: ((^1 + 1, . . . , rife + 1), {wi, . . . , Wk)) 

a 1-homological model of the CDGA j^Ai, Ai being the 

exterior algebra E{xi,ni), if Ui is odd and P{xi,ni) if is even. 



wi = 0, 

for i = 2 to k 

Wi = ds{xi), where the algebra generator of degree Hi 
endfor 



Naturally, the first components of the vector v must be zero, because they correspond 
to the image of the algebra generators with the lowest degree under p. 

Moreover, a general algorithm for computing the 1-homology of CDGAs can be 
described . Clearly, the homology of the 1-homological model obtained can be computed 
using an algorithm based on the establishment of Smith's normal form of the matrices 
representing the differentials at each degree [111, O . 



The computational cost of constructing the contraction {C)s is high. Let us note 
that both the inclusion and homotopy operators of the contraction C^® give an answer 
in exponential time. In fact, the formula of the differential operator ds produced by the 
homological perturbation machinery is given by: 

ds = f5{l-(t)5 + (t)5(t)5 )g. 



With regard to the previous remarks, a first impression is that obtaining ds generally 
becomes a procedure of exponential nature. 

It is possible to take advantage of ds being a derivation. Indeed, the fact that ds 
is a derivation implies that it is only necessary to know this morphism applied to the 
generators of the model (let us observe that there are as many generators as the cardinal 
of the set of indices / indicated). This is an enormous improvement in the computation 
of the differential on the small model. In spite of this, computing ds on an algebra 
generator is extremely time-consuming. 
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4 Inversion Theory 



In this section, we go further in the simphfication of the computation of the differential 
ds- For clarity, we begin this work considering only two algebras. 

As we have seen before, obtaining ds is an extremely expensive procedure. The 
morphism responsible for this is the homotopy operator, 0, due, essentially, to the shuffles 
that are involved in the formulas of (pg^g, and gg^. We intend to eliminate these shuffles, 
and, with this aim in mind, we define the concept of inversion. 



Definition 4.1 Let A and A' be CDGAs and let us consider a homogeneous element 
[oi (g) a'i|a2 (8) • • • lotn ® from B{A<Si A'). We say that a component ^ (8) from that 
element, is responsible for an inversion, if there exists an index j > i with aj 7^ 9 (where 
6 is the unit of A). In this sense, such an element presents k inversions if there exist k 
components responsible for an inversion. 

We will say that an element from B{A^A') has k inversions, if it is a sum of elements 
which each have, as a minimum, k inversions. 

Let us consider the contraction 

(/b«, 9b^, (I>b^) : B{A (8) A') => B{A) ® B{A') 

described in the previous section. We analyze the behaviour of the component morphisms 
of this contraction with respect to inversions. For this purpose, we do not take into 
account the signs in the formulas referred to. 

• The image of an element with at least one inversion under /s^,, is null. 

• The injection qb^, applied to [ai| • • • |an] (8) [a^l • • • |aj„], produces: 

— a unique term with no inversions (that one which comes from juxtaposition), 

— n terms with one inversion, 

— n — 1 terms with more than one inversion. 

• As for the homotopy operator 0^^,, we can state that the image of a homogenous 
element under gives rise to a sum of elements which, if non null, have at 
least one more inversion than the original one. Let us note that an inversion is 
produced by the component a'^ *a' • • • O'n-qi which is always on the left side of those 
components a^i, . . . , a„_, of each summand in the formula (j)B^. 
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Let us consider the contraction which provides us with a 1-homological model for 
the tensor product of two CDGAs, A and A': 

(/, g, 0) : B{A ® A') B{A) ® B{A') hBA ® hBA' (4) 

where 

/ = (/sA ® fBA')fB^y 
9 = 9l3<s{9l3A® 9l3A') 

(f) = 4>B® + 9b,s{4>BA ® 9BA'fBA' + ^BA ® 4>BA')fB,S 

Let us note that the image of an element with an inversion under / is also null, since 
the first morphism applied is fg^. 

Now we assume that there is a perturbation p of the tensor product of the algebras 
A and A'. This perturbation induces, in a natural way, a perturbation 6 on B{A (g) A'). 
Let us analyze the behaviour of such a morphism with respect to inversions. 

Lemma 4.2 Let us consider a perturbation 6 for B{A^ A') induced by a perturbation- 
derivation p for A® A' such that p{A) C A. The image of a homogeneous element with 
k inversions under 6, is a sum of elements with at least k — 1 inversions. 

Proof. 

Let us point out that a component of a homogeneous element from B{A A') is 
responsible for, at most, one inversion and that 6 acts only on a component of the 
element at each term of the resultant sum. 

□ 

Attending to the Basic Perturbation Lemma, one can obtain from the contraction 
(I), a new contraction: 



if 5, gs, ■■ BiA^'A') {HBA ® HBA', ds) . 
We recall the formula for ds'- 

ds = fS{l-(f)S + (f)S(j)S )g. 
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We can observe that / is the last morphism apphed in the formula. If at any stage, 
an element y obtained by applying 0, has more than one inversion, then 5{y) will have at 
least one inversion. In this way, each time we apply 5 0, we obtain a sum of homogeneous 
elements with at least one inversion, and, therefore, the image of these elements under 
/ is null. This means that we only have to consider the summands of having, at most, 
one inversion. 



In consequence, wc can establish the following theorem where we considerably reduce 
the complexity of the computation of ds- 



Theorem 4.3 The formula for (f) , that is involved in the definition ofds, is the following: 

= 0B® + gB®{(t>BA ® gBA'fsA' + 1 ® (t>BA')fB^, 

where 

• 0s®([ai ® oil ■ ■ ■ (g) a^]) 

V- (-l)'^("'^''^)6K_„+i*^---a„)K<8)a;| I" 



0<p<n— g— l<n— 1 

I a' * 4/ • • • a' I <^fi I • • • I Cin-a W \ ' ' 'W 

I n n—q\ I I 9 1 n— q+1 1 I n. 



|a^_i 09 a-_i 



being n — n — p — q and 

'P(ri,p,q) = n - 1 + \[ai\ ■ ■ ■ \an^i]\t + \[a[\ ■ ■ ■ \a'^_g]\t 

p k 

+ X] X] I'^n-q-kl \a'n-q-l\ ■ 

k=o e=o 

• • • • \an] <8) Kl • • • \a'^]) 

= [ai\ ■ ■ ■ \an\a[\ ■ ■ ■ \a'J 

n-1 

+ ^(_l)lk+i|-MI|all [a^i . . . \a,\a[\ai+,\ ■ ■ ■ |a„|4| • • • \a'J 

i=l 

+(_l)KII[ai|-K]| [a[\ai\---\an\a'^\---\a'J. 
Let us note that now the number of summands in the formula above for 0sig, is 

n—l n—q—l 9 i 



2 ' 
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in contrast to the original number of summands: 




2 



in+l 



n-2. 



On the other hand, the formula for gg^ is reduced to n summands, instead of 



( 



m + n 



n 



) 



This theorem is easy to generalize, by induction, to the general case of a twisted 
tensor product of CDGAs ®fgjy4j with I = {l,...,n}, where Ai is an exterior or a 
polynomial algebra with generator Xi and < Therefore, as we saw in Section 

^, the following semifuU algebra contraction can be established: 



where hBAi is a polynomial or a divided power algebra. 

The key to understanding the generalization is the fact that the inversions in 
are those of the last tensor product (as they were defined at the begin- 
ning of the Section, with A = ^'^Z^Ai and A' = An) along with those of = 
i®7=i^i) ® An-i with respect to the last tensor product, and so on. 

Summing up, we obtain an algorithm Alg2 having the same input and output as the 
Algorithm Algl of the last section, but speeding up the steps concerning the image of 
the algebra generators under ds. 

Procedure 2 Algorithm Alg2. 



Input and Output.- the same as in Algl. 



wi = 0, 

for i = 2 to k 

Wi = ds{xj), where Xi is the algebra generator of degree rii 
( using Theorem 4 ) 
endfor 



5 Complexity 

In this section we give a comparison of the algorithms Algl and Alg2 from the point of 
view of their complexity. We are mainly interested in measuring the efficiency of the cor- 



13 



Table 1: Required time 











in Algl 


in Alg2 


s 


= 3 


r = 


2 


0.1 sec. 


0.002 sec. 


s 


= 3 


r = 


3 


0.5 sec. 


0.009 sec. 


s 


= 4 


r = 


2 


31.29 sec. 


0.04 sec. 


s 


= 5 


r = 


3 


3.19 days 


17.6 sec. 


s 


= 6 


r — 


2 


1.45 years 


1.17 min 


s 


= 6 


r — 


3 


24.75 years 


19.56 min 



responding steps concerning the obtention of the differential ds- We consider the degree 
of the algebra generator as the size of an instance. We take as elementary operations 
those ones generating each homogeneous term produced by the different morphisms and 
a worst-case analysis of the algorithms is carried out. 

We calculate the total number of elementary operations needed for computing rf^ on 
a generator of degree k, for both Algl and Alg2. We hold that this number for Algl 
is 

i=0 \j=l J 

and for Alg2 is 

i=0 \j=l 

where ko = mini<i<„|xi| and r is the maximum number of summands given by p{xi), 
where Xi ranges over all the algebra generators of A. 

In the following table, the required time for computing ds{xk) is showed, supposing 
that our computer carries out 10^ elementary operations per second. Let us denote 
s = [k/ko + 1\. Note that, for example, that s = 5 means that if ko — 10, the degrees of 
the algebra generators range over the set {10, 11, ... , 54, 55}. 




6 Implementation Performance 

The algorithm Alg2 has been implemented. The user supplies an encoding of a finite 
CDGA in the form of a twisted tensor product of exterior and polynomial algebras to the 
program, which computes an encoding of a small 1-homological model of this algebra. 

This program is written in Mathematica 3.0, consisting in 300 lines of code and 10 
basic functions. 
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In order to give some indication of the implementation, we report on the time taken 
to compute 1-homological models for certain CDGAs: 

1. E(xi,X2, X3; 1) ® P(x4, X5; 2) E(x6; 3) ® P(xt, 4) ® P(x8; 6). 
Input: ((1,1,1,2,2,3,4,6), 

III I I I I /> ■ /y» /y» <-y» <-y» <-y» <-y» I <-y» <-y» I /y* nf nf rf rf \ 1 

Output: ((2, 2, 2, 3, 3, 4, 5, 7), (0, 0, 0, x^ - x^, x^, 0, x.x^ - 2x^, 0)) 

Time: ds{x4) in 0.55 sec. 
ds{xQ) in 0.27 sec. 
dsi^e) in 0.33 sec. 
dsix-j) in 2.14 sec. 
ds{x^) in 0.28 sec. 

2. E{xi- 1) P(x2; 2) ® P(x3; 6) ® P(a;4; 10) ® P{xr,- 26). 

Input: ((1, 2, 6, 10, 26), (0, — 2xi, a^iXj, 3a;ia;2a;3, 8x1X2X3X4)). 

Output: ((2, 3, 7, 11, 27), (0, -2x,, -8x?, -192xf , 21799895040xP)) 

Time: ds{x2) in 0.16 sec. 
ds^Xj^) in 0.60 sec. 
ds{xj^ in 2.36 sec. 
^{(xg) in 1571.47 sec. 

3. £;(xi; 1) ® P(x2; 2) P(x3; 4) ® E{x4; 5) ^ ^(xg; 7) (0 P(x6; 14). 

Input: ((1, 2, 4, 5, 7, 14), (0, — Xi, XiXa, 0, 2x1X4, — X1X4X5)). 

Output: ((2, 3, 5, 6, 8, 15), (0, -Xi, -x?, 0, 0)) 

Time: ds{x2) in 0.17 sec. 
dsix^) in 0.33 sec. 
ds{x4) in 0.26 sec. 
0^5 (xs) in 0.17 sec. 
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For each algebra, we summarize the results and the time taken to compute a descrip- 
tion of the model. All CPU times are in seconds and calculations were carried out on a 
Pentium III, 128Mb RAM, 7.2Gb Hard disk space. 

This program produces as output an encoding of a certain differential graded algebra 
which could be introduced into another program in order to calculate the homology of 
such objects. We intend to tackle this task in the near future. 
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